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Abstract 

We evolve virtual photon parton densities up to the SUSY threshold and 
higher using coupled inhomogeneous DGLAP differential equations. Reliable 
input parameterizations were available from the c-quark threshold. Limited 

( target photon virtuality ) dependence is observed. The difference to the 
photon structure function is shown to be significant with the introduction of 
SUSY dependent splitting functions. A negligible difference is observed by 
letting the gluino mass enter after the squark mass. An effort is made to in- 
clude the squark threshold effect in such a way that both the renormalization 
group equations are satisfied and the perturbative calculation is reproduced. 



1 Introduction 



There has recently been a great deal of interest in the structure function 
of the photon. This is obtained from the scattering cross-section between a 
highly virtual photon with large square momentum used as a probe and 
a nearly on-shell target photon with square momentum P^, (P^ ^ Q^)- If 
the target square momentum is also large (whilst maintaining the inequality 
P^ -C Q^), the entire structure function can be calculated using renormaliza- 
tion group improved perturbation theory, whereas for low P^ one is limited 
to a determination of the dependence and, as in the case of deep-inelastic 
electron-proton scattering, one needs input information on the structure func- 
tions at some (low) value of Q^, which cannot be determined by perturbation 
theory. A study of the photon structure function as a function of P^ there- 
fore provides information on the extrapolation between the perturbative and 
non-perturbative regimes of QCD. 

Heuristically, one talks about separating the structure function into "di- 
rect" and "resolved" contributions. The former being exactly calculable in 
perturbation theory and the latter involving the uncalculable probability 
that the photon splits into other fundamental particles before being probed. 
Whereas such a picture is useful at the leading order level, higher order cor- 
rections mix these contributions. A formal and more precise analysis was 
first proposed by Witten who pointed out that in an operator product 
expansion for photon-photon scattering the set of operators used in the case 
of photon-proton scattering must be augmented by a tower of photonic op- 
erators, whose matrix elements with the target photon are of order unity. 
In the (more intuitive) DGLAP approach one argues that since the proba- 
bility of finding a particle other than a photon inside a photon is of order 
aemi the probability of finding a "photon inside a photon" is unity plus cor- 
rections of order aem- The DGLAP analysis must then be augmented by 
further off-diagonal splitting functions Kq and Kg which are interpreted as 
the perturbative probability for a photon to emit quark or gluon with a given 
fraction of its longitudinal momentum. 

Interest in photon structure functions has recently been rekindled by the 
prospect of a future high-energy electron-positron collider with centre-of- 
mass energy of up to 1 TeV. Such a machine would enable an investigation 
of the photon structure function over a sufficiently wide range of and 
P^ to provide a stringent test of the evolution of these structure functions. 
Moreover, if the postulated existence of Supersymmetry (SUSY) turns out 
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to be vindicated, these structure functions will reflect the existence of su- 
persymmetric partners within photons. The contribution to the structure 
function due to the crossing of the threshold for the production of squarks, 
was first calculated by Reya 0. However, a consistent analysis of the effect 
of supersymmetry on the photon structure function requires a full analysis 
of the enlarged DGLAP formalism in which above the SUSY threshold the 
standard splitting functions are augmented with splitting functions involving 
squarks and gluinos. This paper reports on such an analysis and displays re- 
sults in which it can be seen that SUSY gives rise to a measurable increase in 
the evolution of the structure photon structure function above threshold. 
Care must be taken to ensure a consistent treatment of the threshold be- 
haviour for squark production as one passes through the threshold and this 
is discussed in detail. 

The outline of the paper is as follows: In Section 2 we discuss the for- 
malism, outlining the extension of the evolution equations to the regime in 
which squarks and gluinos are excited. We also give a description of the 
threshold treatment. Section 3 displays our results obtained from numerical 
solution of the enlarged evolution equations. We show the dependences on 
the SUSY threshold and on the square momentum of the target photon 
as well as on the usual variables and Bjorken-x. In Section 5 we discuss 
our conclusions. 

2 Formalism 

We follow the formalism of Gliick and Reya 0. We will initially be concerned 
with quark and gluon distributions up to the SUSY threshold. 

The nonsinglet (T) and singlet (S) sectors are treated separately. 
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where u, d, s, c, b, t refer to the relevant quark distributions. The factor of 2 
accounts for the anti-quark distribution since for a photon Qi = Qi- f runs 
over aU relevant quark flavours. Each quark distribution is zero at and below 
its threshold hence each new non-singlet (T) is equal to the singlet (S) at 
threshold. 

The evolution is controlled by the following inhomogeneous DGLAP dif- 
ferential equations, 



dT 



Ptt ®T + Kt 



dlnQ'^ 

for each singlet (T) and the coupled equations, 
dE 



dlng2 



-Pee (8) S -|- P^g 



(2.1) 



dG 



GS 



GG 



G + K, 



G 



dlnQ^ 

for the singlet (S) and gluon (G) sector. 

For each distribution F{x, Q^) above, the convolution 



where Pij{x, Q^) consists of the splitting functions p\f in Leading Order (LO) 
and plf in next to leading order (NLO), 



(2.2) 

is defined as, 

(2.3) 



P = 



a.. 



27r 



(0) , 
Pij + 



27r 



(2.4) 



The main difference between the evolution of the photon structure func- 
tion and that of the proton structure function is the presence of the inho- 
mogeneous Ki terms in the evolution equations. Essentially these consist of 
7 quark and 7 gluon splitting functions. They appear in the evolu- 
tion equations without any convolution with a parton distribution since the 
"photon density" inside a photon is taken to be 5(1 — x) up to corrections of 
order ap™,. 



Oie 



27r 



a. 



27r 



27r 



(2.5) 



3 



F2 in LO is given by, 

= {^Ns + (e^) S} (2.6) 

where 

/ J f 

OisiQ"^) evolves according to 

asm ^ 1 _ A In(lngVA^) 
47r /?olnQVA2 /?o' (lngVA2)2 ^ ' ' 

where (3^ = 11 — 2//3 and jdi = 102 — 38//3. All expressions refer to the 
MS renormalization scheme hence we use A;^ which depends /. We evolve 
in NLO to the t-quark threshold and then we evolve in LO thereafter. This 
is because we can only evolve in LO above the SUSY threshold so in order 
to compare F2 for different values of the squark mass Mg we must evolve 
in the same way from the t-quark threshold. Quark masses are taken as 
M(c) = 1.5 GeV, M{h) = A.bGeV, M{t) = 17AGeV. 

Our input data were parameterizations of virtual photon parton den- 
sities taken at the c-quark threshold. A c-quark mass of 1.5 GeV limits 
to less than 1.8 GeV which gives a small ratio r = P^/Q^ ~ 10~^ at high 
Q^. We could not find reliable parameterizations at higher that were 
dependent on P^. 

In order to evolve to the SUSY threshold we use LO and NLO splitting 
functions [^] in (Eq. |2.4| ) and inhomogeneous terms in (Eq. |2.5| ). 

Above the SUSY threshold Ms we are also concerned with squark and 
gluino distributions. As before we have a nonsinglet (S) and singlet (F) 
sector, 



53 = 4.{Sa-Sd) 

S8 = 4{Su + Sd ~ 2Ss) 

S15 = 4{Su + Sd + Ss ~ 3Sc) 

524 = 4{Su + Sd + Ss + Sc - ASb) 
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S35 = A{Su + Sd + Ss + Sc + Sb-5St) 

f 

i 

where Sa, Sa, Ss, Sc, Sb, St refer to the squark distributions. The factor of 4 
arises because Sqf- — Sq^ — Sqf — Sqf. For simphcity all squark distributions 
start at zero at the SUSY threshold M^, although wc could introduce them 
in steps as before with the quarks. The gluino distribution starts at zero at 
the gluino threshold Mg. 

The evolution is controlled by the following inhomogeneous DGLAP dif- 
ferential equations. Each set of nonsinglets are coupled ie. T3 with S3, T8 
with S8, etc... 

dT 

: (g) T + Pts ®S + Kt 



d InQ 

dS 



lnQ2 - PsT^T + Pss^S + Ks (2.8) 

Given that in general the gluino mass Mg is greater than the squark mass 
Ms the nonsinglet sector evolution is given in the region Mg >Q^> M^ by, 



d ing2 

dG 
dlnQ"^ 

dV 



dlnQ2 

and in the region Q'^ > M^ by, 
dE 



Pes S + Peg «) G + Psr ® r + 



= Pgs S + Pgg (8) G + Per r + 



E + Pro (8) G + Prr r + (2.9) 



dlnQ^ 



Pes E + Peg ® G + Psr r + Pel L + Kj: 



dG 
dlnQ2 



= Pgs ® E + Pgg (S) G + Per (^T + Pgl ® L + Kq 
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dr 



d lnQ2 



dL 



d lnQ2 



Prs ® S + Pre ® G + Prr ® r + PrL ® ^ + 



Plj:®T. + Plg®G + Plt®T + Pll®L + Kl (2.10) 



where L is the gluino distribution. In analogy with the Ki used below the 
SUSY threshold, Ky and are the splitting functions of a photon to a 
squark and gluino respectively. In the limit where the gluino mass Mg is 
taken to be equal to the squark mass Mg we do not need (Eqs. p.9|) . 

The Pij{x, Q^) are now a different SUSY set of LO splitting functions 
However in order to reproduce the squark threshold condition, 



Q'' '^ ' > (2.11) 

X 

we use the SUSY set if this applies or the standard set if it does not. Impor- 
tantly this now means that the threshold depends both on x and Q^. This is 
possible because a convolution (Eq. p.3|) only feeds each distribution in the 
region y > x. Hence for a particular Q^, SUSY splitting functions are used 
below a certain value of x and standard splitting functions are used above 
this X during each convolution. 

Our choice of LO inhomogeneous terms (Eq. p.5| ) is made in order to treat 
the squark threshold in a meaningful way. 

The tree level squark contribution to F2 (Appendix A) is important in 
determining the 7 — * squark splitting function i.e., the squark inhomogeneous 
term in (Eqs. |2.8| , ^.9| , p^.lO ). The choice is dependent on the squark threshold 



condition (Eq. p.ll| ), which is function of both x and Q^. At a particular 
there will be a region x > Xs where squarks cannot be produced. 

The contribution to F2 from squark production ( Eq. |A.1| ) is obtained 
from ordinary perturbation theory and includes a term proportional to ln((5^/4M|), 
which is already accounted for as K^, by the use of the full SUSY set of 
splitting functions discussed above. In order to reproduce the correct renor- 
malization group improved evolution far above the threshold we subtract this 
term (apart from terms proportional to / which are negligible far above 
the threshold) from ( Eq. [A.l| ) and introduce the contribution as SB^ in ( 
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Above the SUSY threshold F2 is then obtained from 
\f:J = {gNS + (e^)s} + {S^NS + (e^) 5S }+ 2x3/(e^) 

(the factor of 2 accounts for right- and left-handed squarks). 

The SB^ term is (up to an overall constant) the contribution to F2 from 
squark production, with the removal of the above-mentioned term propor- 
tional to ln(Q7^s)- i-e- 

« F^';,,^,,,- 2x(l-x) ln(QV4Mf), (2.13) 

We note that the difference between using (Eq. [A.l|) and (Eq. [A.2| ) is 
negligible in our case because we are limited to < 1.8 GeV^ at the c- 
quark threshold, giving an r ~ 10~® above the SUSY threshold. 

However this is a different way of treating the threshold behaviour from 
that in 0. At ^ 4M^ it satisfies the Renormalization Group equa- 
tions since the dominant part is in the inhomogeneous term. In the region 
q2 ^ this approach will reproduce the perturbative calculation with the 
correct threshold behaviour up to (ln(Q^/4M^))^. There should of course be 
a small mismatch at large and large x. However we have eradicated this 
by coding the threshold behaviour into each convolution as a choice of SUSY 
or non-SUSY splitting functions. The problem left is that of the discontinu- 
ity at the threshold x boundary. An equivalent problem is discussed in 
involving the DIS^ factorization scheme. Perturbative instabilities mean we 
would have to go to NLO and eliminate the SB^ term by altering what would 
be our NLO inhomogeneous terms (Eq. |2]^). However NLO SUSY splitting 
functions are beyond the scope of the present work. 

Finally it should be noted how quickly F2 changes away from the thresh- 
old with decreasing x. In (Eq. |A.1| ), the term 

v= [l-AM^^x/Q\l-x)Yl^ 

moves rapidly away from zero in decreasing x meaning that the coefficients 
of 

In f — — \ and v 

VI - vj 

in (Eq. [A.l| ) give rise to a real threshold contained in the term SB^ in the 
perturbatively stable x region. 



Up 



SB^ (2.12) 
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To summarize, we take parameterizations of quark and gluon distributions 
inside a virtual photon at the c-threshold. Using DGLAP inhomogeneous 
differential equations we evolve the relevant non-singlet, singlet and gluon 
distributions up to the SUSY threshold. From here we run the distributions 
separately, including or not, the effects of squarks and gluinos. At some \fQ^ 
we form for the virtual photon in such a way as to take account of the 
SUSY threshold condition. 

3 Results 

The variable parameters of the evolution are the (target virtuality), 
(squark mass), Mg (gluino mass), (incident virtuality) and Bjorken x. We 
took these in the ranges, 

< < 1300 Mer 

175 Gel/ < M, < 300 GeV^ 
2m GeV < Mg < ?,mGeV 

700 GeV < ^ < 1500Gel^ 
and in all cases F2 /aem is actually plotted. 

Figure |^ shows a generalised evolution to 1000 GeV . 

The base reference is at lib GeV. There is a considerable difference with 
the SUSY mixing. Note that allowing the gluino mass to be greater than the 
squark mass produces a negligible effect. Note also that the evolution graphs 
coincide above the squark threshold condition (Eq. |2.11|) , this being due to 
it being encoded into each convolution. 

The discontinuity at the threshold can be seen. This results in a sud- 
den drop in F2 below the control non-SUSY graph, this being due to the 
perturbative instabilities already discussed. 

It is important to note the rapidity with which F2 increases as x de- 
creases away from the threshold. This effect due to the SB^ term in the 
perturbatively stable region is discussed in section 2. 

Figure ^ shows dependence up to 1300 MeV . From here on we plot 
Mg = Ms since we have shown the Mg > Ms difference to be negligible. 
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We are limited by our parameterizations in that they are restricted in 
at the c-quark threshold. However non-negligible differences can be noted in 
low X even at = 1300 MeV. 

Figure |^ shows Mg dependence between 175 GeV and 300 GeV. 

The base graph is without the SUSY contributions. As expected, as the 
squark mass Mg increases the distribution's approach to the non-SUSY limit. 
Also the thresholds move to lower x as the threshold condition (Eq. |2.11|) 
requires higher to produce squarks of higher mass. 

Figures ^ and |^ show how F2 varies with a/Q^ at two fixed values of x. All 
graphs show how the distributions must approach the non-SUSY distribution 
for high Mg. However for x = 0.66 we can see the gradual approach to 
a threshold in increasing Mg. For Mg = 275 GeV it is evident that for 
low squarks cannot be produced and the distribution coincides with the 
non-SUSY distribution. Then apart from the perturbative instability the 
distribution rises in higher Q^. 
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Comparative Evolution for ¥ = 1.3 GeV, = 175 GeV, = 175, 300 GeV, Vq' = 1000 GeV 



C\J 
LL 



SUSY Threshold 
Without SUSY 
With SUSY, Mg = 175 GeV 
SUSY, M =300 GeV 




0.2 



0.4 



0.6 



0.8 



Bjorken x 



Figure 1: Comparative Evolution of Structure Function with and without 
SUSY sphtting functions. Difference due to a higher gluino mass Mg is 
neghgible. 



10 



LL 



Vp' Dependence at = 175 GeV, Vq' = 1000 GeV 



Vp' = OGeV 
Vp' = 0.65 GeV 
Vp'=1.3GeV 
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0.4 0.6 

Bjorken x 



0.8 



Figure 2: v dependence of structure function for fixed squark mass Mg at 
a fixed probe virtuality ^/Q^ 
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Figure 3: Dependence of structure function on squark mass Mg at a fixed 
target virtuality and probe virtuality ^/Q^ 
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Q Dependence for fixed X = 0.33 



1.4 
1.3 
1.2 
1.1 
1 

0.9 
0.8 
0.7 



Without SUSY 
M3=175 GeV 
M3 = 225 GeV 
M„ = 275 GeV 



0.6 



600 



800 



1000 



1200 



1400 



1600 



Vq' ( GeV ) 



Figure 4: a; = 0.33 
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1.8 



Q Dependence for fixed X = 0.66 
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Figure 5: x — 0.66 
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4 Conclusions 



We see from Figs. 1-5 that if one can build a machine for which values of 

approach 1 TeV^ (about twice the squark production threshold) there 
is a substantial increase in the value of F2 for the photon. Indeed, the 
evolution between the SUSY threshold and 1 TeV is more than doubled if 
SUSY particles, taken to have a mass of 175 GeV, are present. The difference 
between the structure functions with and without SUSY in the middle range 
of Bjorken-a; is over 30%, which should be easily detectable. 

The effect at = 1 TeV^ is, of course, diminished if the SUSY threshold 
is increased. However, we note that taking the squark masses to 300 GeV 
only has a small effect on ¥2- Conversely, if the squark masses turn out to 
be substantially lighter than 175 GeVQ , (which is not currently ruled out), 
there would be a significant effect on the structure functions at values of 
significantly below 1 TeV^. 

The effect also diminishes if the target photon is off-shell, as will usually 
be the case. However, we see from Figure 2 that this effect is modest. 

The results are fairly insensitive to the mass of the gluino. This is not 
surprising as the gluino contributes very indirectly - it can only be produced 
by a secondary emission from the target photon and then only probed through 
a further interaction with squarks. Taking the mass of the gluino below that 
of the squark, would have had a negligible effect as it is clear that it is the 
squark threshold and not the gluino threshold that must be crossed before 
there is any effect on the photon structure function. 

In summary, we see that the effect of SUSY on the photon structure 
function provides a further good reason for designing a large electron-positron 
collider that would be capable of reaching values of above the SUSY 
threshold for the middle range of Bjorken-x. 

Appendix A 

The contribution to F2 of a left or right handed squark in deep inelastic 
scattering on a photon has been calculated 

^The lowest value we take for the squark mass is 175 GeV since this is above the 
threshold for t-quark and we find it useful to make comparisons of the evolution of the 
structure function in the presence of SUSY with that without SUSY but with six active 
flavours. 
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1-v 



1 + v 



) 



(A.l) 



where, 



r 



V 



[1 - TX^ 



/{i-xr^ 



We calculated this expression for the case ^ 0, where r = P^/Q^. The 
above relation is recovered on r — > with. 



+ B{M^/Q^){l/FG){-48x^r^ + 48x^r + 4xV - Sx^) 

+ B{1 /FG){-12x^r^ + 12a; V - 2x\) 

+ fi(l - 6xV + _ 

+ ln(F/G)(M,Vg2)2(5/r^)(8a;V6) 

+ ln(F/G') (Mf/Q^) (5/r/) (24a;VV& + 2a;V/6+ 12x2 - 2x/b - 2x) 

+ ln(F/G)(B/r;)(i + 6a;V/6+ 12xV/6- 12a;V/6 

-12x^r/b + llx^r/b - 3xV + Ax^ /b - Gx^ + -xr/b 




2 




(A.2) 



where. 



b = 



1 — 2a;r 



F = 



1 + 77(1 - 2a;r) 
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G = 



1 — ri{l — 2xr) 



AM^x 
Q'^{1 — X — xr) 



T] = — Vl — 4:x'^r 
b 

These equations are important in determining the 7 — ^ squark sphtting 
function i.e., the squark inhomogeneous Ks and Kr ■ Also in determining the 
extra SB^ term in (Eq. |2.12|) . The choice is dependent on the squark threshold 
condition (Eq. which is function of both x and Q^. At a particular 

there will be a region x > Xt where squarks cannot be produced. 
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